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Measurements of meson mixing and CP violation parameters obtained by the CDF and DO exper- 
iments at the Fermilab Tevatron are presented. These include results on B s and D meson mixing, 
and searches for CP violation in the decay B + — > J/ipK + , in mixing through semileptonic B s meson 
decays, and in the interference between mixing and decay in the process B s — > J/ifxj)- 
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The Tevatron is a heavy flavor factory in which charm and beauty hadrons of all types are produced at high rates. 
| Even though the hadron collider environment increases the complexity of events, both the CDF and DO experiments 
are making critical measurements of meson mixing and CP violation parameters. CDF exploits its large level one 
trigger bandwidth and level two impact parameter trigger to collect large samples of long-lived particles, and uses 
dE/dx and time-of- flight measurements to suppress backgrounds. DO benefits from its large acceptance for muons 
O ' to collect large semileptonic samples. In this paper, a number of unique results are presented: the measurements 
O ■ of B s and D meson mixing parameters, and searches for CP violation in the decay B + — > J/i/jK + , in mixing 
through semileptonic B s meson decays, and in the interference between mixing and decay in the process B s — > J/tjxf). 



Throughout this paper, charge conjugated modes are implicitly assumed. 



II. MESON MIXING 
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Because the neutral meson mass eigenstates can be a linear combination of the weak eigenstates, they can mix. 

Taking for example B mesons, we can write [l| 
GO 

\B LM )=p\B°}+ q \B° f ), (1) 



where Bjj,l are the mass eigenstates and S2,B9 are the weak (or flavor) eigenstates. The mass eigenstates then 



have masses mn,L — M ± 4p and lifetimes Th,l = T± y. Because of the mass difference, a state produced as 
B°(B°) can oscillate between B° and B° and decay as B°(B°), with a characteristic oscillation frequency Am. In 

1 addition to the variables Am and Ar, it is usual to define 

_ _ i 

Am Ar 

* = —> y=zr (2) 

and 

x = x cos 5 + y sin S, y = —x sin S + y cos 6, (3) 
which are useful when treating integral oscillation probabilities. Here 6 is a strong phase. 



A. B s Meson Mixing 

From the point of view of Feynman diagrams, B meson mixing is dominated by box diagrams with W bosons and 
up-type quarks in the internal lines. Since Vtb is large, diagrams with top quarks are dominant and the oscillation 
frequency is proportional to m\ jm^ . The high oscillation frequency of B mesons was in fact the first experimental 




FIG. 1: Reconstructed B s candidate mass in a hadronic (left) and semi-leptonic (right) channel at CDF. 



indication of a very large top quark mass [2, lia] , and this process is in general sensitive to loop- level contributions 
from new heavy particles. Since V ts > Vtd, B s meson oscillations are faster than Bd oscillations and the frequency is 
correspondingly more difficult to measure. 

The analysis strategy for the measurement of the B s meson mixing frequency, and most of the other measurements 
described here, consists in two main parts. On the "reconstructed side", we have the B s meson under study: its 
flavor at decay time is given by the charges of the decay products, and it is then sufficient to measure its momentum 
and transverse decay length to determine its proper lifetime. To tag the quark flavor at production time, either the 
"opposite side" 6-hadron is used, exploiting the fact that at the Tevatron the vast majority of b quarks are produced 
in pairs, or the "same side" hadron from fragmentation is identified. In the case of B a meson production, this is 
likely to be a K meson. 

The B s meson reconstruction proceeds through the reconstruction of the decay chain B s — > D S X, where X can be 
leptons, or one or three pions. Multiple D s decay channels are considered. While the semi-leptonic channels allow 
for a higher trigger efficiency, they suffer from higher backgrounds and a significantly worse B s meson momentum 
determination due to the escaping neutrino. This is seen in Fig. [T] where the reconstructed B s candidate mass is 
shown for a hadronic and semi-leptonic channel used in CDF. 

The flavor of the B s meson at production is determined using opposite-side tagging variables that include lepton 
charge, jet charge and secondary vertex charge, and same-side tagging based on the charge of a nearby hadron. Both 
experiments use multivariate techniques to combine the various opposite- and same-side tagging variables: CDF uses 
a neural network and DO uses a likelihood. To verify the performance (efficiency and dilution) of the flavor tagging 
algorithm, the Bd oscillation frequency is measured 0, 0| and compared to the world average. 

To measure the oscillation frequency, it is necessary to extract all the information on an event-by-event basis, 
giving events weights in the extraction of the result that are proportional to their probability of originating from a 
certain source. In both experiments, an unbinned maximum likelihood fit is used. Figure [2] shows the measured B s 
meson oscillation frequency in the form of a scan over the likelihood function output. CDF has significantly better 
sensitivity thanks to much larger level one trigger bandwidth, a secondary vertex trigger at level two, and particle 
identification based on dE/dx in the tracker and a time-of-ffight detector. This allows CDF to exploit the excellent 
resolution of the hadronic channels; they measure 0] Am s — 17.77± O.lO(stai) ± 0.07 (systyps^ 1 with a signal larger 
than 5<7. In DO, where the signal is dominated by the semileptonic channels, the signal has a significance of 3cr and 
the measured value is [7J, |8| Am s — 18.53 ± 0.93(stat) ± 0.30(syst)ps _1 . 
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FIG. 2: Measurement of the B s meson mixing frequency. The result is shown in the form of a scan over the value of the 
likelihood function for CDF (left) and DO (right). 
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FIG. 3: Left: dominant "long range" contribution to D meson mixing. The process goes through an intermediate virtual it or 
K meson. Right: doubly cabibbo-suppressed decay. 



B. D Meson Mixing 

As opposed to B meson oscillations, D mesons are expected to oscillate very slowly in the standard model. This 
is because the box diagram, with AC = 2, is very small due to a combination of CKM and GIM suppressions. The 
dominant contribution comes from the "long range" (2 x AC = 1) diagram shown in Fig. [3l and is too small to make 
the oscillation observable before decay. The effect of mixing can however possibly be seen through its interference in 
the doubly Cabibbo-suppressed D° — > K + n~ decay (right panel of Fig. [3]). 

The analysis 0] looks for a mixing-induced time dependence of the ratio of (D° — » K + tt~) / (D° — > K~n + ) decays 
in a sample of events collected using CDF's impact parameter trigger. In this case, the flavor of the D° meson 
at production is tagged by selecting the decay D* + —* ir + D°,D° — > Kir chain. To reduce backgrounds, particle 
identification (PID) is applied as in the B a mixing analysis, and Kir candidates that fall into the D° mass window 
with both mass assignments to the reconstructed tracks are rejected (PID separation has limited discriminating 
power). Events are then selected using the Am = m(KTnr) — m(Kir) — mfir) distributions shown in Fig. [U 

Figure[5]shows the results of the analysis. On the left, the ratio of doubly Cabibbo-suppressed over Cabibbo-favored 
decays as a function of proper decay time exhibits the parabolic dependence 

R(t) = R D + VR^y't + x ' 2 + v ' 2 f (4) 

expected in the presence of mixing, where x' , y' are the mixing parameters defined in section|nl The x 2 /dof for the 
parabolic (flat) fit is 19.2/17 (36.8/19) The right plot shows the probability intervals in the {x l2 ,y') plane. The best 
fit is 3.8cr away from the null hypothesis. 
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FIG. 4: Am = m(Knn) — m(Kn) — m(7r) distributions for Cabibbo-favored (left) and doubly Cabibbo-suppressed candidates 
(right). The combinatorial background and signal are shown in light and dark grey respectively. 




FIG. 5: Ratio of doubly Cabibbo-suppressed over Cabibbo-favored decays as a function of proper decay time (left), and 
(bayesian) probability intervals for the mixing parameters (x' 2 ,y') (right). 



III. CP VIOLATION 

In this paper, CP violation measurements are split in three categories: 

• CP violation in the decay amplitudes, where T(B — > /) ^ T(CP(B) — > CP(f)). In charged meson decays, 
this is the only possible manifestation of CP violation, and a measurement in the B + — > J/ipK + channel is 
described. 

• CP violation in mixing, which can be detected through an asymmetry in charged-current neutral meson decays. 
Recent results from inclusive same-sign dimuon events and exclusive B s — > n + vD~ X decays are given. 

• CP violation in the interference between mixing and decay can be tested in decays with and without mixing 
to a single final state. Here, studies in the B s — > J/ip<t> channel are presented. 




FIG. 6: Left: B candidate mass distribution showing the data and contributions extracted from a fit to the data. The shapes 
of the various contributions are extracted from the simulation. Right: kaon reconstruction asymmetry as a function of kaon 
momentum. 



A. CP Violation in Decays 

DO selects a B + — > J/ipK + (TT + ) sample where the J/tp decays to a pair of muons. This yields a very clean 
sample with approximately 40,000 candidates as shown in the left panel of Fig. [B] To determine the asymmetry, the 
sample is divided in 8 subsamples based on solenoid polarity (3, sign 7 of the J/ipK + system pseudo-rapidity 77, and 
kaon charge q: 

n 0i = ijVe^l + qA)(l + qiA fb )(l + jA det )(l + qpjA ro ){l + q<3A qf3 )(l + faA^), (5) 

where N is the total number of events in the sample, e 13 is the fraction of integrated luminosity taken with solenoid 
polarity (3, A is the CP-violating asymmetry, A fb is the forward-backward asymmetry and all other asymmetries are 
potentially generated by the detector. The total number of events in the J/tpK + and J/ipir + channels, the fraction 
of integrated luminosity for each solenoid polarity and the asymmetries are determined from a fit using the number of 
events in each subsample for inputs. The forward-backward symmetry, as well as all detector asymmetries are found 
to be compatible with zero. At this stage, the observed CP asymmetry needs to be corrected for the asymmetry in 
kaon reconstruction efficiency induced by the different interactions of K + and K~ with the detector, which is made 
of matter only. This reconstruction efficiency asymmetry is measured from a sample of D* + — > D°ir + , D° — > fivK~ 
events as a function of kaon momentum (right panel of Fig. [S]), and applied to the kaon spectrum in the B + sample. 
After this correction, the CP asymmetries are measured to be Acp(B + — > J/'ipK + ) = 0.0075 ± 0.0061(siat) ± 
0.0027(syst) and A C p(B+ -> J/ipn + ) = -0.09 ± 0.08(stat) ± 0.03(syst). The dominant systematic uncertainty 
comes from the mass distribution model. 



B. CP Violation in Mixing 

CP violation due to mixing between neutral B mesons can be searched for by measuring the asymmetry in same-sign 
dimuon events: 

A w = N{bb -> fi+fi+X) - N(bb -> f£j£X) 
SL N(bb^ ii+n+X) + N(bb^ fJ.-fJ.-X)' {> 

Two approaches are used: 



• CDF llj uses all dimuon events but exploits the two-dimensional (/i — fi) impact parameter significance dis- 
tributions to unfold contributions from different sources, including beauty and charm hadrons. CDF finds 
A sl = -0-0080 ± 0.0090(siai) ± 0.0068(syst). 
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FIG. 7: B s —> J/^4> candidate samples after event selection in DO (left) and CDF (right). 



• DO [12j uses all dimuon events, and estimates the contributions from all possible processes to the sample, includ- 
ing sequential decays, Drell-Yan, instrumental, etc. The "8 subsamples" technique described in section IIII Al 
is used to extract the asymmetry. DO measures A^ L = —0.0053 ± 0.0025(sioi) ± 0. 0018 (syst). 

In both cases, the asymmetry from fake muons due to kaons and pions which have asymmetric reconstruction 
efficiencies (see section UlI A[) is determined from data and corrected for. 

The flavor-specific asymmetry A S SL then needs to be derived using the sample composition and the known oscillation 
probabilities, since Ag^ has contributions from both Bd and B s mixing. This is done by taking A d SL from the B- 
factories, and the known Bd and B s production ratios and mixing parameters. The results are: 

• CDF: A S SL = 0.020 ± 0.021(stat) ± O.OW(syst) ± 0.009 (inputs); 

• DO: A S SL = —0.0064 ± 0.0101(a^Z uncertainties combined). 

An alternative approach is to use flavor-specific decays, as is done in the DO analysis [3| of B s — > n + vD~X, where 
again the asymmetry is determined using the 8 subsamples technique. This yields A S SL — 0.0245 ± 0.0193(s£ai) ± 
0.0035(syst). 



C. CP Violation in the Interference Between Mixing and Decay 



The search for CP violation in the interference of mixing and decay can be searched for in final states that are 
common to both B s and B s decays. Both CDF 14 1 and DO [3] make flavor-tagged measurements of the CP 
properties of the B S (B S ) — * J/ip4> process. This is similar to the measurement of sin(2/3) in Bd(Bd) — * J/tpK s decays 
at the B- factories. The angle 4> s (= — 2/3 s ) is however much smaller. An easy way to see this is that the opposing 
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, whereas for (3 it is proportional to 
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side in the corresponding unitarity triangle is proportional to 
This measurement requires the determination of the CP eigenvalue of the vector-vector J/ip4> final state, and benefits 
from the identification of the flavor of the decaying B s meson. 

Both experiments use events where the J /ip meson decays to a pair of muons, leading to a high trigger and 
reconstruction efficiency. At the event selection level, CDF uses a neural network with variables that include PID 
from dE/dx in the drift chamber and information from the time-of-flight detector. DO uses a "square cuts" event 
selection without PID. The resulting B s meson candidate samples are shown in Fig. [3 
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FIG. 8: Two-dimensional confidence intervals in the (/3 S ,AT S ) plane from the flavor-tagged analysis of B a — » J/ip(f> decays 
from DO (left) and CDF (right). 



Since the J/ifxj) final state consists of two vector mesons, both CP eigenvalues are possible. The value for a given 
B s candidate is measured by determining the relative polarization of the J/tp and 4> mesons: the angular dependence 
of the relative directions of the decay products fi + and K + K~ is expressed in the J/ip restframe ( "transversity 
basis") in terms of this polarization. The corresponding probability density functions are determined from simulation, 
and used to determine each candidate's probability to be in a given eigenstate. Note that the detector efficiencies 
are not flat as a function of transversity angles, and this is properly accounted for. 

Constraining the flavor of the decaying B s meson removes a twofold ambiguity in the result. This is achieved by 
determining the B s meson flavor at production and using the mixing frequency measured by CDF (see section THAI) . 
Both experiments use both opposite-side and same-side tags as described in section lll Al The tagging power obtained 
in CDF is eD 2 = 1.28% and 3.65% for opposite- and same-side tags respectively, and 4.68% for the combined tag in 
DO. 

To maximize sensitivity, an unbinned likelihood fit is used which assigns event-by-event signal and background 
probabilities based on reconstruction uncertainties, dilution, accuracy of CP eigenstate determination etc. DO uses 
constraints on strong phases from world-average values [3] of measurements of B^ — > J/ipK* decays, whereas CDF 
lets these phases float. Constraining the phases essentially eliminates the remaining ambiguity in the extraction of 
the result, but there may be differences between the phases in Bd — > J/ipK* and B s — > J/ipcf) decays. From the 
likelihood fit, the average lifetime, Ar s , 4> s (= — 2/3 s ), the magnitudes of the polarization amplitudes and the strong 
phases are extracted. 

A small added complication to the extraction of the final CDF result is due to the fact that in that analysis the 
likelihood profiles are not parabolic close to the minima. Since standard frequentist techniques can't be applied, CDF 
uses a Feldman-Cousins-like likelihood ratio ordering to build the two-dimensional confidence region in the ((3 S ,AT S ) 
plane. The results obtained by CDF and DO are shown in Fig. [8] 

Both experiments measure values that are in reasonable agreement with the standard model: the "standard model 
probability" is 15% and 7% for the CDF and DO results respectively. However, it is clear from Fig.[8]that both results 
pull in the same direction (4> s (= — 2/3 s )) and with similar magnitudes. While an official combination of these results by 
the experiments is not ready yet, the UT/i£ collaboration has performed a combined fit of these results with the CP 
asymmetries described in section llllBI They try multiple methods to unfold the strong phases constraint from the DO 
result, and in each of the methods they find a value of 4> s that is 3<r or more away from the standard model prediction. 
It should be noted that for both CDF and DO, the measurement uncertainties are completely dominated by statistics, 
so that the results will still improve by quite a bit, hopefully leading to a conclusive statement. If the UT fit conclusion 
is confirmed (and strengthened), this would be the first evidence for CP violation outside the CKM mechanism . 



IV. CONCLUSIONS 



The beauty and charm physics programs at CDF and DO continue to produce a large number of excellent results. 
A number of recent highlights were presented here. CDF has now seen evidence for D meson mixing, and both 
CDF and DO are making precise measurements of CP violation asymmetries in B meson mixing and decay. The 
B a meson mixing parameters are being pinned down: CDF measures Am s with high precision, both experiments 
find a value of Ar s in good agreement with the standard model, and both experiments see similar deviations of <p s 
from the standard model expectation. The latter measurement will become significantly more precise with increasing 
statistics. If the result is confirmed, this will be the first sign of CP violation in the quark sector outside the CKM 
mechanism. 
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